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We report progress in our exploration of the finite-temperature phase structure of two-flavour lat- 
tice QCD with twisted-mass Wilson fermions and a tree-level Symanzik-improved gauge action 
for a temporal lattice size A^^ = 8. Extending our investigations to a wider region of parameter 
space we gain a global view of the rich phase structure. We identify the finite temperature tran- 
sition/crossover for a non-vanishing twisted-mass parameter in the neighbourhood of the zero- 
temperature critical line at sufficiently high j3. Our findings are consistent with Creutz's conjec- 
ture of a conical shape of the finite temperature transition surface. Comparing with NLO lattice 
;i;PT we achieve an improved understanding of this shape. 
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1. Introduction 

The goal of the tmfT collaboration is to explore the applicability of the 2-flavour Wilson- 
twisted-mass (Wtm) fermion formulation set-up as described in Ref. for investigations of lat- 
tice thermodynamics. We refer to [^] for the definition of the Wtm action and its parameters 
K and jJ.. The staggered-fermion formulation is computationally inexpensive however con- 
ceptually controversial [Q]. On the other hand, the non-perturbatively improved Wilson fermion 
formulation requires the calculation of the improvement coefficients and furthermore needs oper- 
ator improvement. Thus, the Wtm formulation, nowadays combined with a tree-level Symanzik- 
improved gauge action [^, appears as an appealing alternative for finite temperature lattice sim- 
ulations, whose potential should be explored. It offers automatic 0'{a)— improvement by tuning 
the bare, untwisted quark mass only. This makes high-statistics simulations for thermodynamical 
problems affordable with pseudoscalar masses as low as m^i > 300MeV [|6p. 

As a necessary preparatory step we characterise the phase structure of the model by locating 
the transition/crossover lines and surfaces in the j8 — k: — /x— coupling space. In contrast to our 
report at Lattice 2007 we now find clear evidence for the Aoki phase at j8 < 3, lower than 
searched for previously, and a first order phase transition surface in an intermediate jS region be- 
ginning with j3 f» 3.4, consistent with earlier findings. For larger couplings 3.6 < jS < 3.8, we 
observe a thermal transition at Kt{T ^ 0) near Kc{T = 0), which emanates from the first order like 
structure at smaller j3 -values. With growing j3, it splits up to become a succession of two tran- 
sitions confinement^deconfinement— ^confinement in the K"-direction. The confinement aspect of 
these transitions is exhibited by the Polyakov loop and its susceptibility, the chiral aspect by the 
pion norm (the integrated pseudoscalar correlator) and the scalar and pseudoscalar condensates. 

Our ultimate goal is to investigate the finite temperature transition (or crossover) at maximal 
twist and physically relevant pion masses. The most obvious procedure of directly scanning the 
transition at maximal twist has proven the hardest to tackle, since its precise location is yet un- 
known. Still, our efforts so far gave us more detailed insights than reported in into the phase 
structure at non-vanishing twist parameter /i for a larger range of j8 -values. Our findings are qual- 
itatively consistent with predictions from ;^PT and continuum symmetry arguments 

2. Preview of the phase structure 

There are three sources for our expectations concerning the phase structure of 2-fiavour Wtm 
LQCD at finite temperature: (i) The phase structure at zero temperature was predicted using 



;^PT [Q, |10|, |11|] and subsequently verified numerically [|12|]; the parameter space at T = was 



shown to contain a phase of broken parity-flavour symmetry (Aoki-scenario) in the j3 — K"— plane 



at strong coupling [ |13[ ] (for Wilson gauge action) and a first order phase transition surface in the 
j8 — /I— plane (the first order scenario of Sharpe and Singleton) adjacent to the former and extending 
into the weak coupling region. Both scenarios crucially depend on the choice of fermion and gauge 



action [14]. (ii) The general properties of the phase structure for (non-) perturbatively improved 



Wilson fermions, as established by e.g. the CP-PACS [15] and DIK [16] collaborations, should be 
found in the j8 — K"— plane at weaker coupling, too. (iii) Based on the argument that in the continuum 
theory the physical transition temperature cannot depend on the twist angle, Creutz [j8|] conjectured 
the finite temperature phase transition/crossover surface to be closed and form a cone (cf. the right 
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panel of Fig. |I|) around the critical line K"c(j8) for a range of j8 -values somewhere between the Aoki 
phase and the deconfinement transition in the quenched limit (jSdec ~ 4.5 at K" = 0). The location 
of the tip of this cone and the detailed shape of its surface are subject to lattice artifacts, as may be 
the nature of the transition for different twist angles. 
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Figure 1: Left: Map of our simulations performed at /i < 0.007. Right: Emerging phase structure; regions 
identified by tmfT measurements are marked by magenta symbols. 

3. Simulation results 

The simulations were mainly performed on lattices of spatial size A^fj = 16 and Nt = 8. The 



algorithm in use is described in Ref. [17|. An overview of the regions covered by our simulation 
is presented in the left panel of Fig. [l[ It contains all transition and crossover points that we found 
in our simulations and is supplemented by data describing the T = critical line fQ;(j8 > 3.75) 
provided by the ETM collaboration. 

After the unsuccessful search reported in we have checked the existence of the Aoki phase 
at the three values j3 G {1-8, 3.0, 3.4}. The necessary and sufficient condition for spontaneous 
symmetry breaking reads 

lim lim (v//75tV>a'<„/i / 0- (3-1) 



Figure S shows the Fisher plots [ 1 3 ] of the order parameter ( y) ^ function ofh/{ xjfiYs t y) > 
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Figure 2: Fisher plots for the order parameter {xj/iysT^ xj/) describing its hmit for h~ 2kh^0; left: j3 = 1.8, 
right: (5=3.0. 

with /j = 2k"/i, for j8 = 1.8 and 3.0. The condensate has a finite limit for ^ (a positive in- 
tercept at /j = 0) for suitable K values. For the smallest value j3 = 1.8, the left panel of Fig. ^ 
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gives convincing evidence for the onset of spontaneous symmetry breaking somewhere in the range 
0.24 > K" > 0.23. The right panel strongly indicates that at j3 = 3.0 the symmetry-broken phase - 
if it exists at all - is realized only in the interval 0.203 < K" < 0.204. A volume extrapolation would 
be welcome but is unavailable at this j8— value. This marks the lower boundary K"]°*'^''(j8) of the 
Aoki phase. We have not systematically searched for its upper boundary K'"^'"''^(j8). 
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Figure 3: Signals of metastability at j3 = 3.4 and 3.45. The two-state signal left: in the pseudoscalar 
condensate, right: in the real part of the Polyakov loop. 

At j3 = 3.4 we reinforce our preliminary results published in Ref. |0] that, instead of an Aoki 
phase, there are strong metastabilities when simulating at finite twist /i « 0.0068. This confirms 
the first-order Sharpe-Singleton scenario, and represents a remnant of the phase structure at T = 
(on symmetric lattices). This is shown for j3 = 3.4 and 3.45 in Fig. |3|for the order parameter (left) 
and the Polyakov loop (right). The behaviour of the latter at this transition cannot be interpreted as 
thermal. Also the average plaquette (as checked at j8 = 3.4 and K = 0.1827) is running in separate, 
metastable histories starting from hotter/colder neighbouring states. 

In the range 3.4 < j8 < 3.65 between strong and weak coupling, simulating aljx = 0.0068, we 
see a finite temperature transition, located at Kj ^ 'Q(j8), i.e. far above the chiral line for T = 0. 
The transition is inherited from the first doubler structure existing at higher K. The corresponding 
Kt{P,IJ.) evolves rapidly with j3, moving closer towards the critical line JQ.(j8) with increasing jS 
(cf. Fig. [I]). The transition is possibly of first order as indicated by histograms of the Polyakov loop 
at j8 = 3.6 in the interval 0.18 < K" < 0.24: they show a relatively flat but clear double peak structure 
at K" = 0.22. On 32^ x 8 lattices simulated at the same jS — K" — /i— point one observes long-living 
metastability in the Polyakov loop. Also at higher j3 = 3.75 this thermal transition continues to 
exist, still separated from the unfolding thermal transition surface around (j3). Finally, it appears 
to join the latter in the neighbourhood of j8 = 4.0. Before this joining happens. Fig. ^ shows the 
characteristic behaviour of the Polyakov loop and its susceptibility measured at jS = 3.75 with 
jj. = 0.005 over an interval 0.16 < K" < 0.21 that is covering both transition regions. The two peaks 
are clearly visible, but the lower susceptibility peak is not yet well resolved on the coarse k scale 
of Fig. I 

It is our goal to work at maximal twist, thus ultimately we need to resolve the lower peak in 
Fig. ^. The position and nature of the first (i.e. lower in k) finite temperature transition/crossover, 
when studied at 7^ 0, changes considerably between the metastability region described by Fig. || 
and the scaling region. Midway between the two, at j3 = 3.65, one sees the Polyakov loop entering 
and then leaving again a narrow deconfining region in K (cf. the left panel of Fig. B|). Approaching 
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Figure 4: Left: real part of the Polyakov loop, right: its susceptibility as functions of k for j3 =3.75 and 
II = 0.005. The dashed line marks Kc{P = 3.75, T = 0). 



the physically relevant scaling region, p > 3.75, the peak structure grows even clearer. The narrow 
peak (supposed to contain the deconfined phase) moves gradually further to lower K with rising 
j3 = 3.75, 3.775 and jS = 3.8 (cf. the right panel of Fig. |). 
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Figure 5: The real part of the Polyakov loop as function of k; left: for j3 =3.65 with ji = 0.0068, right: 
for J3 = 3.75, 3.775 and 3.8 (all with n = 0.005). 

With a large computational effort the rise and decline of the Polyakov loop begins to be re- 
flected also in the corresponding susceptibility. The Polyakov loop susceptibility gradually de- 
velops a double peak which is already clearly visible at jS = 3.8 (cf. the left panel of Fig. |6|). 
Remarkably, the lower— k peak is accompanied by a peak of the pion norm (presented in the right 
panel of Fig. ^ emphasizing the entrance into the deconfining and chiral-symmetry restored phase. 
So far, within the presently available statistics, we could not detect a significant similar chiral signal 
marking the exit. 

4. The cone signal and comparison with lattice 

We interpret the structured signal in the Polyakov loop near outlined in the preceding section 
as a numerical indication for passing through a confinement^deconfinement crossover approach- 
ing K"c from small positive quark mass, followed by a deconfinement^confinement crossover at 
small negative quark mass. This picture is qualitatively consistent with the conical structure pre- 
dicted by Creutz (see also Fig. |I|). The thermal transition takes place at a given value of the quark 
mass which is determined by the parameters K and jJ.. At tree level the according relation reads: 

2 



1 



1 



(4.1) 
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Figure 6: Left: the susceptibiUty of the Polyakov loop, right: the pion norm as functions of K for )3 = 3.75, 
3.775 and 3.8 (all at ^ = 0.005). 



Using lattice chiral perturbation theory (L;^PT), one obtains the following modification of (4.1) in 
NLO (cf. [0]): 



Z2' 



Z|4 



K K 



(4.2) 



From the zero temperature simulations of the ETM collaboration we know that fQ.(j3 = 3.75) = 
0.1660(1), Zs « 0.6 and Zp « 0.3. K is an unknown 0{a) coefficient that is introduced by L;^PT. 
The twist angle (O is defined by tan O) = /i / (0.5 ( 1 / — 1 / jq) ) . 

In Fig. ^ the thermal transition points that we have found at j8 = 3.75 for jj. = 0.005 and 
0.007 are compared with the two relations above. In principle, K and niq can be determined from 
a fit. However, given the rather large uncertainties of Zp and Zs and also the small number of 
available data points, we are only able to check whether the formulae are capable of describing the 
data. For the tree level formula this is not the case (in the plot: niq = 0.01, JQ- = 0.1660), due to the 
K (2/Kc — 1/ K")-symmetry of equation (11). However, choosing = 0.028, K = 0.5, Zs = 0.6, 



Zp = 0.3 and jq. = 0.1660 the NLO formula (4.2) is consistent with the data points 




Figure 7: Thermal transition points found at j3 = 3.75 with — 0.1660 compared to tree level (m^ = 0.01, 
green) and LxVT (m^ = 0.028, K = 0.5, Zs = 0.6, Zp = 0.3, blue) curves. 



5. Summary and outlook 

We investigated the global phase structure of 2-flavour Wilson-twisted-mass LQCD with a 
tree-level Symanzik-improved gauge action. Our simulation results support the existence of an 
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Aoki phase (at p < 3.0) as well as the first-order Sharpe-Singleton-scenario at finite twist (for 
j3 3.4). Moreover, we find indications of a conical structure of the finite-temperature crossover 
surface in a bounded j3— interval, 3.65 < j8 < 3.8, with /i = 0.005 and 0.007. Our next steps will 
be the localization of additional finite temperature crossover points to refine the L;^PT— prediction 
of the transition line and to bound the value /ij- for the thermal crossover at maximal twist within 
a narrow interval. Depending on the findings of this investigation, a subsequent simulation at 
maximal twist to detect /I7 numerically will follow. 
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